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We derive the single molecule equation of state (force-extension relation) for DNA molecules 
bearing sliding loops and deflection defects. Analytical results are obtained in the large force limit 
by employing an analogy with instantons in quantum mechanical tunneling problems. The results 
reveal a remarkable feature of sliding loops - an apparent strong reduction of the persistence length. 
We generalize these results to several other experimentally interesting situations ranging from rigid 
DNA-protein loops to the problem of anchoring deflections in AFM stretching of semiflexible poly- 
mers. Expressions relating the force-extension measurements to the underlying loop/boundary de- 
flection geometry are provided and applied to the case of the GalR-loop complex. The theoretical 
predictions are complemented and quantitatively confirmed by MD simulations. 
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In nature DNA is rarely found in its straight "naked" 
state as usually depicted in the introductory pages of 
elementary textbooks. In most "in vivo" situations an 
overwhelming fraction of DNA is rather strongly configu- 
rationally constrained by binding proteins causing loops, 
bends and deflections. The advent of single molecule 
stretching techniques Q has opened the possibility of 
measuring the " equation of state" of single tethered DNA 
molecules in a variety of different conditions 0- While 
the statistical mechanics of unconstrained DNA under 
tension is theoretically well understood in the framework 
of the Worm Like Chain model Q the presence of topo- 
logical constraints like supercoiling 0, 01 or geometrical 
constraints like protein induced kinks and bends (Qj, |3 
renders analytical results more difficult. In this letter we 
expand the repertoire of analytically solvable " equations 
of state" by deriving the force-extension relation for a 
DNA molecule featuring loops and large deflections in 
the limit of strong stretching forces (for the small forces 
case see 0). The computation is performed by evaluat- 
ing quadratic fluctuations around the looped solution - 
a non-constant saddle-point of the DNA elastic energy. 
The method is essentially analogous to the semiclassical 
treatment of tunneling amplitudes in quantum mechan- 
ics and instantons in quantum field theory 'lC)|. After 
deriving the general result (that we accompany with MD 
simulations) we focus on two interesting experimental ap- 
plications: the stretching of the GalR loop complex 
and of tangentially anchored semiflexible polymers from 
a surface in AFM experiments. 

Stretching a sliding loop. In the following we neglect 
the DNA twist degree of freedom for it is not constrained 
from outside, and no external torsional torques are act- 
ing on it. In this case DNA of length L is described 
in the continuum limit by specifying only the unit vec- 
tor tangent t(s) to the chain. Here s is the contour 
length along the DNA with -L/2 < s < L/2. The 
chain is submitted to an external constant force so that 
the kinetic plus potential energy of such a chain reads 
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FIG. 1: Examples for large deflections in a stretched DNA 
molecules, a) A freely sliding linker protein stabilizes a DNA 
loop, b) A rigid ligand with opening angle a causes a kink 
in the DNA. c) Tangentially anchored DNA stretched from a 
surface by an AFM tip. The tilting angle r as well as anchor- 
ing angles a and a can strongly affect the elastic response. 
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ds with the bending stiff- 



ness A = IpksT where Ip is the orientational persistence 
length and ksT the thermal energy; for DNA at room 
temperature Zp w 50 nm [l2| . 

In the following we parametrize the tangent as t = 
(cos0cosi?, sin0cosz5,sin-(5) and put the force along 
the X-axis so that the potential energy part writes 
— F cos (j> cos -d . Note that the angle is measured with 
respect to the equatorial plane (as on a globe). This 
parametrization is necessary to take properly into ac- 
count the inextensibility of the chain imposed by the con- 
dition = 1. In the following it is convenient to intro- 
duce the dimensionless contour length t = s/X with the 
deflection length A = ^ AjF 0, • The latter becomes 
the relevant lengthscale characterizing the loss of orien- 
tational correlation in the case of DNA under large ten- 
sion (replacing the usual tension-free persistence length 
Ip — A/kBT). In these coordinates the chain energy 
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The relevant saddle point - the loop in the x-y plane 
shown in Fig. la - satisfies 5E{) = with the boundary 
conditions "dioop (±L/2A) = 0, (±L/2A) = and 2tt. 
In the limit of large lengths i/A — > oo the loop solution is 
the well known "kink" solution from theory of quantum 
tunneling 10] cos (pioop = 1 — 2cosh~^(t) and i9/oop = 
with the loop energy given by Eioop — 8\/FA — LF. 
For large values of V AF we can expand Eq in terms of 
quadratic fluctuations {S'd, d(j)) around the looped saddle 
point {-dioopAioop) 
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with the in- and out-of-plane fluctuation operators 
t|| = -d^/dt^ + (1 - 2cosh~2 1) and = -d^/dt^ + 
(l — 6cosh~'^i) respectively. A closer inspection of the 
discrete spectrum of the two operators reveals the physics 
behind. Ty has a zero eigenvalue resulting from the 
translational shift-invariance of the loop along the chain 
that costs no energy for L/X oo. The absence of neg- 
ative eigenvalues is in agreement with intuition as in 2D 
the loop is a (topologically) stable saddle point. The out- 
of-plane fluctuation operator Tj^ shows a richer behavior. 
Again it possesses a zero mode, this time resulting from 
rotational invariance of the problem. More remarkably, 
in contrast to T|| the out-of-plane fluctuation operator 

T_L has a negative eigenvalue —3. The latter underlines 
the intrinsic instability of the loop in 3D as described by 
the elastic energy expressions ^ and |2 This brings us to 
the question how to describe the obviously stable phys- 
ical situation shown in Fig. la. In order to model the 
action of a sliding linker we introduce a term account- 
ing for the DNA self-interaction. In lowest order it can 
be written as a function of the perpendicular distance 
Azc = A f_]jm5'd{t)dt w A f_^^J§{t)dt of the two 
overcrossing DNA arms at the contact point tc ~ 1.915 
(resulting from the crossing condition t(, = 2tanhic). 
The total energy of the chain can now be written as 
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with a short-ranged but otherwise arbitrary attractive in- 
teraction potential V acting at the crossing point. Note 
that if V (Azc) has a minimum at Azc = the loop sad- 
dle point stays unaffected by the self-interaction. DNA 
is known to effectively attract itself in many solvents de- 
spite its strong negative bare charge. Typical situations 
inducing DNA self-attraction are poor solvents (like alco- 
hol, small neutral polymers like PEG) 0], the presence 
of multivalent counterions (like CoHex and Spermidine) 



or small cationic proteins acting as linkers between two 
DNA surfaces. Single molecule stretching experiments on 
DNA condensed with multivalent counterions performed 
by several groups jl^ rnight bear loops or related struc- 
tures like DNA toroids [iq. 

The partition function of the system in 
Fig. la can generally be written as Qioop — 
Jioop ^ (i^ - 1) [i] e"'^^'°*[iJ where Ji^^^ represents 
the path integral over the functional neighborhood 
of the loop solution and the (5-function enforces the 
chain inextensibility. This partition function is nothing 
but the Euclidean path integral of a quantum particle 
moving on an unit sphere under the influence of an 
external constant force. Using the 0, parametri- 
sation the partition function can be rewritten as 
Qi^^p = /X>[0]X'[^]e-'5^*°*[''''^l-i"'^" with a metric 
term Cm = exp (J 5 (0) ds log (cos (i^))) resulting from 
the inextensibility constraint. It can be shown that this 
term does not contribute at the quadratic level of the 
approximation because we expand around iJjoop = 0. By 
virtue of the decomposition of the in- and out-of-plane 
fluctuations at the quadratic level (Eqs. |21 and |3Jl the 
partition function can be conveniently factorized 
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and with the in- and out-of-plane plane partition func- 
tions 
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The notation J* reminds that the (translational and rota- 
tional) zero modes of Ty and T_l respectively have to be 
handled with care. A naive Gaussian integration would 
lead to a formal divergence in both cases. After dealing 
with this problem by introducing collective coordinates, 
a method well known from tunneling theory [loj . and 
applying the Gelfand-Yaglom method for the computa- 
tion of the fluctuation determinant we obtain the in-plane 
partition function Q|| = ^PLFe"^^/^ . 

The computation of Q± which follows similar lines of 
reasoning requires rewriting the (5-function in 13 in its 
Fourier representation. After introducing He, the char- 
acteristic function of the interval [—tc,tc] (i.e. Ilc{t) — 
1 for t G [—tc,tc] and = otherwise), we rewrite 
Si!}dt = He (t) (5t? (t) dt as a scalar product. The 
path-integral [7| becomes a Gaussian with a source term 
ifcllc and can be solved by constructing the Green's 
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This leads to the result 



8^^ (x)''' e-^e^^l^"' with r = 2/(9^3 - 30t,) « 0.35 
a numeric constant. Combining those results we obtain 
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The resulting force extension relation (Ax) 
ksT-^ logQioop writes in leading order 
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In the limiting case of a very deep potential V (z) strongly 
localized around z ~ the last term is negligible and 
the force-extension relation becomes independent of the 
detailed nature of the contact interaction. If we in ad- 
dition consider large forces, the O (l//3iF)-term can be 
neglected and the relation can be cast in a more illu- 
minating form 
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that resembles the well known loop-free WLC response 
(Az) /L « 1 - (2y^VF^ 0, but with a strongly 
renormalized apparent "persistence length" l"^^ : 
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Equations EH and show that one has to be very cau- 
tious when interpreting experimental stretching data in 
terms of persistence length and stiffness. The presence of 
a loop modifies the elastic response of the chain in such a 
manner that the persistence length can appear effectively 
reduced as stated in^] For a single loop this is obviously 
a finite size effect involving the scaled total length L/lp. 
But the effect remains significant over a large range of pa- 
rameters: for Ip/L — 0.1 one finds l"^^ w 0.31^p whereas 
for Ip/L = 0.02 there is still a remarkable effect, namely 
0.74/p. 

To check the prediction of Eas.llOlan dllll we performed 
molecular dynamics (MD) simulations '17]. As seen from 
Fig. 2a the theoretical predictions quantitatively agree 
with the simulation results. 

Applications. In order to generalize the asymptotic 
(large F) relations 1101 and [TTl to other interesting situa- 
tions we can rewrite them in a more illuminating fashion 
in terms of the three relevant lengths L, Ip and A 
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FIG. 2: a) Stretching DNA looped with a sliding linker: 
Comparison of MD simulation results (datapoints) with the- 
oretical prediction (solid lines) for the force-extension rela- 
tion as given by Eq. IIUI The force is measured in units of 
Fp — ksT/lp. The different slopes correspond to different 
values of Ip/L. The dashed line shows free chain (no loop) 
behavior, b) Stretching of a GalR loop complex: MD simula- 
tion (dots) vs. theory, Eq. 1131 (solid lines) for various values 
of Ip/L. 



The first term in the brackets is the length-loss by ther- 
mal fluctuations. Because of the small size of the loop in 
the large F limit (A <C L) it comes as no surprise that the 
term X/2lp coincides with the thermal fluctuation contri- 
bution of a loop-free chain. On the other hand the second 
term in the bracket AX/L describes the loss of length due 
to the (loop induced) elastic deflection. This asymptotic 
decomposition of the two contributions in 1 121 leads to im- 
mediate generalizations. After short inspection it is easy 
to see that any strongly localized 0| DNA deflection 
(like in Fig. lb and Ic) can be appropriately continued 
and mapped piecewise onto fractions of the full loop so- 
lution, Fig. la. By this reasoning a deflection angle of any 
size (between and 27r) can be immediately related to 
the length loss in complete analogy to the loop case. 

As a first example consider the stretching of rigid 
DNA-protein complexes which come as large kinks or 
fixed angle loops, cf. Fig. lb. A prominent example is 
the GalR repressosor complex that was studied in single 
molecule experiments |llj . Following the upper reason- 
ing one obtains the same force-extension relation as in 
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1101 but with an apparent persistence length given as a 
function of the opening angle a of the complex: 
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For the GalR complex the opening angle is not known but 
there are some indications for the "anti-parallel" config- 
uration, i.e., a « 

[llII3- WithUniandlTSlat hand one 
can predict the loss of length due to the presence of a 
GalR complex with the conjectured angle a = 0. For the 
force F — 0.88 pN applied in the experiment |0 with a 
loopsize of 38 nm [ll|. Eg s. 1101 and 1131 predict a loss of 
length of 56 nm (with the DNA hidden inside the loop 
included). Remarkably, the experimental value by Lia 
et al. (llj is 55 nm ±5 nm. The latter result together 
with ^] and ^] gives convincing evidence for the " anti- 
parallel" loop model. An independent check ofl 101 and 1 131 
is provided by MD simulations (shown in Fig. 2b) which 
were performed for various ratios of Ip/L. In conclusion, 
Eq. ^1 enables one to directly measure the angle a - an 
important feature of the DNA-protein complex geometry. 

A second application concerns AFM stretching exper- 
iments of semiflexible polymers A stable anchor- 
ing can be achieved when the polymer is tangentially at- 
tached at its two ends as in Fig Ic. Force extension data 
in such a setup have to be interpreted with care. The 
boundary anchoring angles at the AFM tip as well as at 
the surface - a and a respectively - can significantly alter 
the measured apparent persistence length, a fact which 
has been overlooked before. A more trivial effect of a tilt- 
ing angle r in Fig. Ic (between the line of contact-points 
and the force direction) can alter the result in addition 
[2]|. A simple calculation for large forces F again gives 
the same functional relation as in 1101 with an "apparent 
persistence length" given by 
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While the angle r can be completely canceled by shifting 
the tip in the lateral direction, the influence of a and a 
cannot be fully eliminated by this procedure. This means 
that for short to intermediate sized semiflexible polymers 
{Ip/L ~ 1 — 1/30) the influence of boundary conditions 
has to be taken into account through Eq. 

In conclusion, the force-extension relation of looped 
DNA which we derived in the limit of strong stretch- 
ing force can be generalized to DNA featuring large de- 
flections. In particular our approach provides an an- 
alytical expression for the force-extension response of 
DNA-bending proteins which can give important infor- 
mation with regard to the DNA-protein complex geome- 
try. Large deflections due to the anchoring at the AFM 
tip and at the substrate also affect the measured persis- 
tence length for not too long molecules. Finally these 



effects might be related to the strong reduction of the 
persistence length found in condensed DNA stretching 
experiments by Baumann et. al. 15]. 
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